MKPH-T-01-10 

Construction of Consistent Meson Exchange Currents by Laplace Transform^ 



Hartmuth Arenhovel and Michael Schwamb 
Institut fur Kernphysik, Johannes Gutenberg-Universitat, D-55099 Mainz, Germany 

(February 8, 2008) 

We propose a new method for the construction of a consistent meson exchange current in r- 
space for the spin-isospin dependent central and tensor part of phenomenological nucleon-nucleon 
potentials by using a Laplace transformation, which allows the representation by a finite number 
of Yukawa functions. This method is applied to the Paris and the recent Argonne Vig potentials. 
Results are presented for electrodisintegration of the deuteron near threshold. 
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The contribution of meson exchange currents (MEC) to electromagnetic reactions (e.m.) on nuclei, like photo 
induced reactions and electron scattering, constitute an important manifestation of meson degrees of freedom medi- 
ating the strong interactions between nucleons in nuclei. They describe the e.m. interaction with a nucleus during 
the interaction of nucleons and appear formally in the nuclear current operator as two- or many-body operators. 
Thus they are intimately linked to the TV TV-interaction. However, for a given TVTV-potential there exists no a priori 
way of constructing the appropriate MEC-operators, unless the TVTV-potential is derived from an underlying more 
■ fundamental fieldtheoretical framework with explicit subnuclear degrees of freedom which allows the construction of 
ly-^ ] the corresponding two-body currents in parallel to the TVTV-potential. Although the existence of exchange currents 
associated with an exchange TV TV-potential has been acknowledged for a long time, the early realistic potentials being 
phenomenological to a large extent prevented thus the explicit consideration of such exchange currents. 

A breakthrough came with the meson-theory based construction of MEC-operators. Thus purely meson exchange 
models of the TV TV- interaction like the Bonn potentials [Q allow one to build the appropriate consistent MEC-operators 
uniquely (for a recent derivation including leading order relativistic contributions see Q). In these studies it turned 
out that the most important MEC contributions came from the tt- and p-MEC which are directly related to the spin- 
"V* \ isospin dependent central and tensor part of the TV TV-potential ||-|6| . But for largely phenomenological potentials the 
construction of a consistent MEC remained questionable. In view of the fact that realistic phenomenological potentials 
^ | incorporated 7r-exchange as longest range contribution, one often used for such potentials only a regularized 7T-MEC 
. as an approximation. 

However, the above mentioned results for meson-theoretical models suggested a method which provides also for the 
spin-isospin dependent central and tensor part of a phenomenological potential a consistent MEC whose construction 
is based on an analogy with the properties of tt- and p-MEC Essentially this method relies on a splitting of 

the spin-isospin dependent central and tensor part of a given TV TV-potential into a tt- and p-like part for which the 
consistent MECs are known. While the approach of Riska || is based on the momentum space representation of the 
potential and can be applied to any given phenomenological potential, the method of Buchmann et al. [7]] is conceived 
for r-space calculations and needs a representation of the potential as superposition of Yukawa functions. Thus the 
application of the latter method seems to be more limited because not all phenomenological potentials have such 
a form. On the other hand, in view of recent high precision, though largely phenomenological TVTV potentials, like 
the Argonne Vis [|l0| > h would be desirable to have a method available which allows the construction of a consistent 
exchange current for such potentials directly in r-space representation. 

It is the aim of the present brief note to show, that it is possible to represent a given realistic r-space potential in 
general as a superposition of Yukawa functions such that the construction of a consistent MEC in r-space is easily 
achieved. This new method is based on the representation of the radial central and tensor parts of a spin-isospin 
dependent potential by a Laplace transformation. First, we will briefly review in the next section the approaches 
of Sect. Ill contains the central idea introducing the Laplace transformation. Explicit applications to the Paris 
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and Argonne potentials are presented in Sect. [V. In Sect. [V] we consider electrodisintegration of the deuteron 
threshold as a test case for the evaluation of the corresponding MECs. Finally we close with a summary. 



II. BASIC FORMALISM 

We will focus now on the spin-isospin dependent part of a given AW-potential having the general form 

V aT {r) =n-T 2 (ai ■ a 2 V£ T {r) + S 12 ^ r (r)) , 

where 

5i2 = 3 5\ ■ f c?2 • f — <?i • <?2 
denotes the usual spin tensor operator. A pure unregularized 7r-exchange potential is given by 

V: T (r ) = fx • r 2 ai • V a 2 . V J m ^ ( r ) , 
mi 

with the potential strength denoted by V® and 

p -rnr 

J m (r) 



Airr 

By recoupling 

a x ■ V S 2 ■ V = i {ox • a 2 dl c + S 12 d^ T 
with radial differentiation operators 

<c-- r ^(r-) and = r ±{\ ± ■) , 



and using 



with 



®r,C Jrn(r) = m 2 J rn (r) ~S(r) , 
dr T Jm{r) = m 2 F T (mr) J m (r) , 



F T (x) = 1 + - (1 + -] 

X X 



it can be brought into the form (hi) with 



VS(r) = ^dl c J m „(r) 



mi 

= V° F T (m ff r) J m „(r). 
Analogously, the dominant part of the p-exchange potential has also the form (Q) with 

V£(r) = 2 V? dl c J nlp (r) and V£(r) = -V p ° d 2 r , T J rrlp (r) . 

The potential strength is denoted by V® . 

The corresponding exchange currents read for 7r-exchange 



ri,f 2 ) = Tf 2 V° (8(2 - ri) a l a 2 • V 2 J m „(|n - r 2 |) - (1 «-> 2) 

+CTi-ViCT 2 -V 2 J m7r (|ri-x|) V* J m „(\x - ? 2 \)) , (13) 

and for p-exchangc 

j p (x, r l7 r 2 ) = Tf 2 V° (5(x - f x ) (a 2 x V 2 ) x B x J„ lp (| n - r 2 |) - (1 <-> 2) 

+ (cti x Vi) ■ (<? 2 x V 2 ) J mp {\fi -x\) y x J m „(\x - fyj) , (14) 

where we have introduced 

Tl 2 = (n x f 2 ) 3 . (15) 

Defining the currents in momentum space by 

MQUQ2) = 77^ / d 3 xd 3 ri d 3 r 2 e-^ s e^ j(x , f u r 2 ) , (16) 



these currents read 



(2tt) 3 



Jtt {q,quqi) = -iS(q- q\ - q 2 )Tf 2 [at (a 2 ■ q 2 )v n (q 2 ) - (1 <-> 2) 

+ — "4 (^i ' 3l)(^2 ■ 92) (v^(qi) - v n (q 2 ))] , (17) 

qf-qi 

3 P {q,qi,q2) = -iS(q- <fi - g 2 )rL 3 2 [d*l x (<?2 x q 2 )v p {q 2 ) - (1 <-»• 2) 

+ ^ 4 (^i x 9i) ' (^2 x 9*2) K(9i) - ",,(92))] , (18) 



where 



v,/ P {q) = / dW, /p (ry- r (19) 



denotes the corresponding Fourier transforms of the potentials. 

As mentioned in the introduction, for the case of phenomenological central and tensor parts of a realistic spin-isospin 
dependent iV TV-potential, the construction of a consistent MEC is based on the idea to split these potential terms into 
a 7r-like and a p-like part for which the corresponding MEC are known. This method has been developed independently 
by Buchmann et al. 0] and Riska Q for the Paris potential || . Buchmann et al. start from the representation of the 
Paris potential in terms of Yukawa functions which reads in particular for the spin-isospin dependent part 



12 



VS T (r)^Y,9c:,JmAr), (20) 
3=1 
12 

W (r) = E 9%i F t Jm, (r) ■ (21) 



Here rrij, g^J and g^J denote appropriate masses and coupling constants (for details see 0). Using (0) and (|J) one 
finds 



12 CTT 12 (TT 

Vc (r) = %o E §f J ^ (r) + ( E ) ^) ■ (22) 

12 ar 

^(r)=a r 2 T E^ J ™,« (23) 



These expressions can be rewritten as 
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12 err 

VS T (r) = dl c (v*- uk *{r) +2V"- like (r)) + (]T ^)<5(r) , 



m 

3=1 3 



V£ T (r) = dt tT (v*- like (r) - V p - Hke (r)^j 



where the 7T- and p-like parts are given by 

y(7r/p) — like 



12 



with coupling constants defined by 



9i 



3 = 1 

2 gift) and 5 ; = _L( 5 -- 5 -) 



(24) 
(25) 

(26) 
(27) 



3 3 

Thus the representation of the spin-isospin dependent part of the Paris potential in terms of n- and p-like potentials is 
achieved if the (^-function disappears in (|24|), which means the coupling constants have to fulfil the following condition 

12 12 (jt 

9c,j 



3=1 3 = 1 3 



0. 



(28) 



If this condition is not fulfilled, as is the case for the Paris potential, one can modify in (|20j) the short range part of 
Vq t (t) by splitting it into 



VS T (r)=VS T (r)+VS r (r) 



(29) 



such that Vq t (t) is identical to Vq t (t) in the long and medium range part but obeys (p8[ ). In this is achieved by 
changing the coupling constant of the highest mass g c T i 2 — > $c r i2' wnere t ne latter is determined from fl28|), i.e. from 



9c, 12 



i=i 



2 \ ^ 5c j 
12 ^ m? 



(30) 



The difference g c r sr = 5ci2 — 5c~i2 serves as coupling constant for Vg r (r) = g c T sr J mi2 (r). For the remaining short 
range potential Vff" a consistent MEC is easily constructed jj]. It is obvious that this procedure is not unique and 
that it introduces some ambiguity, which indeed is characteristic for such phenomcnological approaches. The hope is 
that the important physics in the long and medium range part is preserved. 

Riska on the other hand considers the momentum space representation of the potential terms obtaining the it- and 
p-like pieces from 



V 



3p 2 



47T 
47T 



where the subtraction of the term 



drr 2 (vS T (r)j (pr) + 2 F^(r)j 2 (pr)) - v°J (0) 
drr 2 (vS T (r)j (pr) - V^(r)j 2 {pr)) - v&(0) 

V ° T (0) = 4tt J drr 2 VS T (r) 



(31) 
(32) 

(33) 



is required in order to eliminate the ^-function. It constitutes again a modification of the short range part of the 
original potential part as mentioned above. In fact, in this case it amounts to the subtraction of a <5-function 

VS r (f) = v c T (0)S(r). (34) 

Also here the above mentioned ambiguity becomes apparent, because any finite range potential V^ r (r) with the 
property 

4tt J drr 2 V£ r (r) = v°7 {0) (35) 

would also serve to eliminate the 5-function. 

The corresponding MECs are obtained from (|l^) and ( |l4| ) by replacing the Fourier transforms v n and v p by 
the corresponding ones, y n ~ hke and v p ~ hke , respectively. This method appears more general than the approach of 
Buchmann et al. since it does not rely on the Yukawa representation. However, we now will show that also the latter 
method can be applied to potentials with a more general radial dependence. 
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III. REPRESENTATION BY A LAPLACE TRANSFORM 



Our first idea was to approximate the given potentials by a series of Yukawa functions whose coefficients and masses 
are obtained by a least square fit. But depending on the accuracy needed this can be quite a formidable task because 
of the high dimensional parameter space involved. But then it turned out that an easier and more systematic approach 
can be based on the representation of the potentials by a Laplace transform. Indeed, because of the fact that realistic 
potentials contain as longest range contribution a 7r-exchange potential, one can represent the central and tensor parts 
of a given realistic potential by a continuous superposition of appropriate Yukawa functions (compare with ( pp[ ) and 

m,i. e . 



Vc(r) = / dmgc(m) J m +m^(r) 
Jo 

poo 

= J m Ar) / dmg c (m)e- mr , (36) 



V T (r) = / dm g T {m) F T ((m + m„)r) J m+m „(r) 
Jo 



= #, r J m » / dm ^ ' e~ mr , (37) 

which is essentially a Laplace transform representation. We note in passing, that the Fourier transform of Vc is simply 
given by 

vc{q) = r dm ^J, ^ V2 . (38) 



Introducing 7r- and p-like potentials as in (24) and fl25| ) with corresponding representations 

/•OO 

V^"- like [r)= dmg„ /p (m)J m+m „(r), (39) 
Jo 

one obtains the following relations 

V c( r ) = f Q dm(g 7r (m) + 2g p {m))({m + m 7 r) 2 J m +mJr)-5(r)y (40) 

/•OO 

Vr(r)= / dm (m + m 7! ) 2 (g 7! (m) - g p (mj) F T ((m + m Tr )r) J m+m „(r) . (41) 
Jo 

Comparison with ( |36] ) and (|37j) gives 

gc(m) = (m + m„) 2 (g w (m) + 2g p (m)) , (42) 
5t(to) = (to + m 7r ) 2 (.g 7r (m) - ,9 p (to)) . (43) 

Again, in order to eliminate the 8- function in (p0|), one needs the condition 

i p oo f \ 

dm{g v (m) + 2g p (m))= / dm , g ^ \, =0. (44) 



If this is not fulfilled, one has to separate again a short range potential as in (29) 

1 f°° 

Vgr(r) = ~ dmg s c r (m)e~ {m+m T )r , (45) 



where the coefficients g^J (to) in principle can be chosen quite arbitrarily except for the fulfilment of the relation 

(m + m^y J (m + m„y 

However, in practice one would choose them such that only the short-range part of the original potential is modified. 
The associated exchange currents read 
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h--iik&(x,n,f 2 ) = Tf 2 J dmg„(m) (s(x - n) 3\ a 2 • V 2 J m +m x (ki - r 2 \) - (1 <-+ 2) 

+(?! ■ Vi (T 2 ■ V 2 J m+m „(|ri - f|) V K J m +m„(|x- raj) J , (47) 

J P -iike{x,ri,r 2 ) = Tf 2 J dmg p (m) (s(x - ft) (3 2 X V 2 ) x 3\ J m+m „(|ri - r 2 |) - (1 <-► 2) 

+ (cti x Vi) • (<7 2 x V 2 ) Jm+m^Qri - x\) Vi J m +m,(|x - r 2 |)) . (48) 



For the explicit application it is useful to discretize the integrals in (|36|) and Q37]) , for example by Gauss quadrature 

AT 

N, (49) 



/'OC J » 



where AT denotes the number of Gauss points, the Gauss points and Wj the corresponding weights. Then one can 
represent the radial functions in ( |40| ) and ( |4l| ) by a finite number of Yukawa functions 

AT 

Vc{r) = w i G^K) + 2 ffp (m,)) ^ 2 (r) , (50) 
3=1 

AT 

VHr) = XI ^ ^(^(m-j) - gpimj)) F T (fijr) J Mj (r) , (51) 
3=1 

where we have set /Zj = mj + m^. The consistent exchange current is then given as superposition of n- and p-like 
currents as listed in ( |l3| ) and ( |T4| ) 

Jrnec = Jn-like{x, ?l, f 2 ) + Jp-like(x, ?1, f 2 ) (52) 

with 

AT 



" 7r _i ifce (f,fi,r 2 ) = rf 2 X g-n{mj)\8{$ - fi)3i 3 2 • V 2 J w (|ri - r 2 |) - (1 <-> 2) 

• Vict 2 • V 2 J^(|fi -x|) Vx J H {\x-r 2 \)), (53) 
(f,ri,f 2 ) = Tf 2 Wj gp(mj) fflpr - fi) (<? 2 x V 2 ) x oi J^(|ri - r 2 |) - (1 2) 

x Vi) • (ct 2 x V 2 ) J^dfi - f|) Vx - f 2 |)) . (54) 



3=1 

AT 

J p— like 

3=1 



In order to determine the unknown coefficients, one can choose an appropriate grid of TV radial points (i = 
1, . . . , AT), and obtains an inhomogeneous set of 2N linear equations for the 2A^ coefficients g v (rnj) and g P (jrij) 



N 



E4(^+ 2 ft)' ( 55 ) 



3=1 
N 



3=1 



where we have introduced for convenience 



Kj = g^rrij), pj = g p (m 3 ),, 

ci = v c (r t )/J m An), u = V T (n)/J m An), (57) 

It is worth mentioning that the matrices Afj and Afj do not depend on the potentials, only on the chosen grids of 
radial points r$ and masses rrij . 



G 



Solving the equations ( p5| ) and ( |56| ) by matrix inversion, one finds for the coefficients 

N 

3 



/°3 

The condition (M4) becomes 



i N 

^E(^ c )7/^-^ T )7/^)- ( 59 ) 



3 * ■■ 



^^■(^+2p,) = 0. (60) 

3=1 

If this is not obeyed for the given central potential, we split off a short range Yukawa potential for the highest mass 
Hn = m N + run 

VS r (r)=c^—. (61) 

This means that the coefficients c, in ( |58|) and (|5^) have to be changed according to 

Ci ^ Ci + ce- mNri , (62) 

leading to new coefficients Ttj and pj. Then the unknown coefficient c is determined by the requirement that the 
resulting new coefficients iij and pj fulfil (pTj|), which leads to 

c = ivT^ T- • ( 63 ) 



IV. APPLICATION TO PARIS AND ARGONNE Vis POTENTIALS 

As first example, we have chosen the Paris potential |J , because in this case we can compare the Laplace transform 
representation of the associated MECs directly with the consistent it- and p-like MECs of . Since for a numerical 
evaluation one needs a reliable representation only for the radial range r = — 10 fm, we have chosen the grid of 
TV radial points Tj in the range between and about 12 fm with variable stepsize with the highest density of points 
close to the origin, where the potentials exhibit the largest variation, and then with increasing stepsize approaching 
the highest value tm- In fact, an educated choice will have to take into account the radial behaviour of the potential 
under consideration. In detail, we have chosen the r-grid to be defined by the following expression 

Tj =r Q + (e a( ^ 1} - 1) e 6 ^- 1 ) fm , (64) 

with ro = 0.01 fm, a = 0.01, and b is determined for a given N by the requirement that the highest point r/v lies 
approximately between 10 and 12 fm. The parameter b is listed in Tab. [|. In order to check the convergence with 
respect to the number of points we have considered N = 12, . . . , 20 in steps of 2. The Gauss points and weights for 
the integral over m have been chosen according to 

7T 

rrij = s tan(— Xj + 1) , (65) 

W > = I cos'd^ + l) ' (66) 

where Xj and yj are Gauss points and weights, respectively, for integration between —1 and 1. For a given N the 
scale factor s is determined by minimizing the mean absolute deviation between the Yukawa representation V^ Y,N ^ 
and the original potential V 

A ( V <y,A0) = 1__ f dr \ V ( Y > N )(r) - V(r)\ , (67) 

Tmax ' J rn 
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evaluated between ro = .01 fm and r max = 10 fm. The resulting scale factors and the relative mean deviations, i.e., 
A(V {Y - N) ) divided by the average potential strength \V\ 

\V\ = 1 — f maX dr\V{r)\, (68) 

for the central and tensor potentials are also listed in Tab. [| 

In Fig. [l] we show the original central and tensor potentials together with their Yukawa representations where 
we have multiplied them by the inverse of the pion Yukawa function J mir (r) for A — 12, . . . , 20 in order to exhibit 
in greater detail the accuracy of the representation. One readily notes the rapid convergence and the very good 
representation over the whole range for A > 16. Indeed, the Laplace transform representations for N = 18 and 20 are 
undistinguishable from the original form on this scale. Therefore, we show in addition in Fig. || for TV = 16, 18, and 
20 the relative deviations of the Laplace transform representation from the original form. For small r-values up to 
about 3 fm the deviations are extremely small whereas for higher r-values the relative deviations become somewhat 
larger, but this is of little importance in view of the rapid fall-off of the potentials themselves with increasing r. It 
is interesting to note that the relative deviations are larger for the central part than for the tensor potential. The 
reason for this feature is the rather rapid variation of the central potential near the origin. 

The analogous results for the Argonne Vis potentials are displayed in Fig. || and Fig. |^ with scale parameters and 
relative mean deviations listed in Table |J. One readily notes again the excellent representation for N = 12, . . . , 20. In 
this case the relative deviations are larger in the tensor part, because the central potential exhibits a much smoother 
behaviour near the origin compared to the Paris potential. In fact, the Laplace transform representations of the 
central part for N = 18 and 20 are undistinguishable from the original form in Fig. ^ and even on the magnifying 
scale of Fig. |] one notes only very tiny deviations. 



V. RESULTS FOR DEUTERON ELECTRODISINTEGRATION NEAR THRESHOLD 



A benchmark for the study of meson exchange current effects in electromagnetic reactions on nuclei is deuteron 
electrodisintegration near threshold at higher momentum transfers P,|ll|-|l4j . The threshold region is dominated by the 
Ml-excitation of the antibound 1 S'o-resonance in ./VA-scattering at very low energies. With increasing momentum 
transfer the inclusive cross section at backward angles, where the transverse current contribution via the inelastic 
transverse form factor dominates, the one-body current contribution drops rapidly due a destructive interference 
of deuteron S- and Z?-wave contributions. In this situation the contribution of MEC, which are of shorter range 
than the one-body currents, becomes relatively more important, in fact dominant. Only inclusion of such MEC 
gives a satisfactory description of experimental data yielding thus clear-cut evidence for the presence of exchange 
currents @ IJ^-R ■ 

We show in FigT^the backward inclusive cross section near threshold at a cm. excitation energy of E np = 1.5 MeV 
as function of the momentum transfer squared, calculated for the consistent ir- and p-like MEC according to Q and for 
the new Laplace transform representation for A = 12, 16 and 20. It is almost impossible to note a difference between 
the different curves on this scale. For this reason we exhibit in Fig. [| the relative deviation between the Laplace 
transform representation for N = 12, 16 and 20 and the orig inal MEC of @. Already for A = 12 the maximum 
deviation does not exceed 1 %, for A = 16 it is less than 0.1 %. In fact, the difference to the A = 20 result is hardly 
noticeable even on the enlarged scale of Fig. Ej. For N = 20 the agreement is perfect. This clearly demonstrates 
that one has achieved already with A = 12 quite a satisfactory parametrization, while for N — 16 an almost perfect 
description for the consistent MEC is obtained. Here, we do not compare to experimental data for which one would 
need to include additional contributions, because we only want to demonstrate that the new method works very well. 

We then have evaluated the analogous MEC contributions for the Argonne Vis potential. Also in this case we found 
an excellent convergence of the Laplace transform representation as is demonstrated in Figs. [?] and || (note the further 
enlarged scale) , where we show the same quantities obtained for the Argonne Vis potential as for the Paris potential. 
One readily notes that in this case the convergence is even more rapid. For A = 12 the maximum deviation from the 
case N = 20 is about 0.2 %, and for N = 16 it is less than 0.01 %. 

We will end this section with a comparison with experimental data. In Fig. |^ the theoretical results obtained for 
the Paris and the Argonne potentials are exhibited together with experimental data from Refs. Jl5| — |lT[| . The theory 
includes besides the consistent tt- and p-like MEC in addition relativistic one-body current and wave function boost 
contributions of leading order in p/M . One readily notes a satisfactory agreement for both potential models with 
experiment up to a squared momentum transfer of about 25 fm -2 . At higher momentum transfers the theory deviates 
significantly from experiment. However, in this region one expects a break down of the present approach in view of 
the applied p/M-expansion [ p"8| . 
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VI. SUMMARY AND OUTLOOK 



In this note we have shown that it is possible to construct directly in r-space a consistent meson exchange current 
for a spin-isospin dependent iV TV-potential by representing the potential as a continuous superposition of Yukawa 
functions, essentially a Laplace transform representation. In this way it is possible to rewrite the potential into a tt- 
and a p-like part, whose corresponding consistent MECs then serve as a basis for a consistent MEC for the given 
potential except for a small modification of the short range part. 

The feasibility of this method by discretizing the continuous superposition into a finite number of Yukawa functions 
has been demonstrated first for the Paris potential for which a consistent r-space MEC exists already. For a given grid 
of N masses the corresponding coefficients of the Yukawa functions are uniquely determined by a properly chosen grid 
of N radial points and involve a simple matrix inversion only. It turned out that the convergence with the number 
of terms is very rapid, and that with N = 16 one obtains an excellent representation of the potential. The same was 
found for the more recent Argonne Vis potential. 

The resulting consistent MEC, represented by a corresponding superposition of tt- and p-like MECs, has then be 
checked by evaluating the inclusive cross section of deuteron electrodisintegration near threshold. For the associated 
observable, the inelastic transverse form factor, which is sensitive to MEC, we found for the Paris potential excellent 
agreement with previous evaluations and for both potentials a very rapid convergence with the number of terms. 
Thus the present method will easily allow one to implement a consistent MEC into an r-space calculation using 
phcnomcnological iV TV-potentials. 

Finally we would like to emphasize that, although consistency of the MEC with the potential is achieved, one has 
to be aware of the fact, that this MEC is by no means unique. We have already alluded to a certain arbitrariness 
in separating a short-range part in order to eliminate an otherwise appearing <5-function. Furthermore, there is in 
addition a freedom in the spin-operator structure as has been pointed out already in [Q . Only the longest range part 
of the 7r-like exchange current, namely the genuine 7T-MEC is on safe grounds. 
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TABLE I. Parameter value b for the radial grid, scale parameter s for the mass grid of Gauss points and relative mean 
deviation for Paris potential as function of the number of Gauss points N. 



N 
b 

s [fm- 1 ] 
A(V^ N) )/\V C \ 
A(V^ n) )/\Vt\ 



12 
0.42 
1.64 
0.25 • 10" 
0.20 ■ 10" 



14 
0.34 
2.20 
0.21 • 10" 
0.35 • 10" 



16 
0.29 
2.67 
0.54 • 10" 
0.72 • 10" 



18 

0.245 

3.09 
0.72 ■ 10~ 7 
0.65 ■ 10~ 6 



20 

0.21 

3.48 
0.14 ■ 10~ 7 
0.12 ■ 10~ 6 



TABLE II. Scale parameter s and relative mean deviation for Argonne Vis potential as function of the number of Gauss 
points N. 

N 12 14 16 18 20 

s [fin" 1 ] 2.19 2.43 2.75 3.04 3.235 

A(V^ Y ' N) )/\V C \ 0.21 -10~ 3 0.69 • 10~ 4 0.58 ■ 10~ 4 0.45 • 10 -4 0.36 • 10~ 4 

A(V^ Y ' N) )/\V T \ 0.61 ■ 10~ 3 0.24 ■ 10~ 3 0.11 ■ 10 -3 0.31 ■ 10 -4 0.98 • 10~ 5 



> 




FIG. 1. Central and tensor spin-isospin dependent parts of Paris potential. Notation: original form: full curves, Laplace 
transform representation for different number N of Gauss points: N = 12: dotted, iV = 14: dashed, N = 16: dash-dot. 
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FIG. 2. Relative deviation of Laplace transform representation of central and tensor spin-isospin dependent parts of Paris 
potential from original form for different number iV of Gauss points: Notation: iV = 16: dash-dot, N = 18: long dashed, 
N = 20: full curves. 
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FIG. 3. Central and tensor spin-isospin dependent parts of Argonne Vis potential: original form and Laplace transform 
representation. Notation as in Fig. [[]. In the central part the various Laplace transform representations are almost undistin- 
guishable from original form. 
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FIG. 4. Relative deviation of Laplace transform representation of central and tensor spin-isospin dependent parts of Argonne 
Vis potential from original form for different number N of Gauss points: Notation as in Fig. bl 
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FIG. 5. Inclusive cross section for deuteron electrodisintegration near threshold for Paris potential for final state excitation 
energy E np = 1.5 MeV and electron scattering angle e = 155°. Notation: solid curve: consistent MEC of Q coinciding with 
consistent MEC of Laplace transform representations for iV = 12, 16 and 20; dashed curve: without MEC. 
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FIG. 6. Relative deviation of inclusive cross sections for deuteron electrodisintegration near threshold for Paris potential 
calculated with MEC from the Laplace transform representation (d 2 a™' ') (N — 12, 16, 20) to the one with MEC from M 



(d a) for final state excitation energy E np 
dashed, N = 16: dash-dot, N = 20: solid. 



= 1.5 MeV and electron scattering angle 9 e = 155°. Notation of curves: N = 12: 
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FIG. 7. Inclusive cross section for deuteron electrodisintegration near threshold for Vi8 potential for final state excitation 
energy E np — 1.5 MeV and electron scattering angle 8 e — 155°. Notation: dashed curve: without MEC; solid curve: consistent 
MEC of Laplace transform representations for iV = 12, 16 and 20, which coincide with each other. 
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FIG. 8. Relative deviation of inclusive cross sections for deuteron electrodisintegration near threshold for Vis potential 
calculated with MEC from the Yukawa representation (d 2 a (Y ' N ' > ) for TV = 12 and 16 to the one with TV" = 20 for final state 
excitation energy E np = 1.5 MeV and electron scattering angle 8 e — 155°. Notation as in Fig. pi 
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FIG. 9. Inclusive cross section for deuteron electrodisintegration near threshold for Paris and Argonne Vis potentials for 
final state excitation energy E np = 1.5 MeV and electron scattering angle 8 e = 155°. Notation: dotted curve: without MEC; 
solid curve: consistent MEC for Argonne Vis potential; dashed curve: consistent MEC for Paris potential; Experiment: filled 
circles: @, open circles: @ (0 e = 155°, averaged over energies OMeV < E„ p < 3 MeV); filled triangles: @ (0 e = 180°, 
averaged over energies OMeV < E np < 10 MeV). 
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